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The direct CP violation parameter Re(e'/e) has been measured from the decay rates of neutral kaons into 
two pions using the NA48 detector at the CERN SPS. The 2001 running period was devoted to collecting 
additional data under varied conditions compared to earlier years (1997-99). The 2001 data yield the result: 
Re(e/e') = (13.7 ± 3.1) x 10~ 4 . Combining this result with that published from the 1997,98 and 99 data, an 
overall value of Re(e'/ e ) = (14.7 ± 2.2) x 10~ 4 is obtained from the NA48 experiment. 



1. INTRODUCTION 

1.1. CP violation in the neutral kaon sys- 
tem 

CP violation has been discovered in the neutral 
kaon system in 1964 |IJ. The main component 
of the effect occurs in the mixing between 
Kq and K . The physical states K$ and Kl 
deviate from pure CP = ± 1 eigenstates, with 
the mixing described by the parameter e. Di- 
rect CP violation can occur in kaon decays to 
two pions through the interference of amplitudes 
with different isospins ||. This is described by 
the parameter e'. The quantity which can be 
measured experimentally is the double ratio R of 
the decay widths: 

T(K L -» A°)/r(ifg -> A ) 

T{K L -> TT+TT-)/T{K S -> 7T+71-) 

« 1 - 6 x Re(e'/e) 

In the Standard Model, CP violation arises 
from the irreducible complex phase in the CKM 
matrix Q. Direct CP violation is predicted by 
the Standard Model, with typical computations 
for Re(e'/e) ranging from » -10 to 30 xl0~ 4 §. 

From the data taken in 1997-98-99, NA48 pub- 
lished § a result (15.3±2.6) x 10~ 4 . KTeV § 
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published a result based on the data taken in 
1996-1997 (20.7±2.8) x 10~ 4 . These recent results 
showed the existence of direct CP violation in the 
neutral kaon system. 

We report here on the measurement of Re(e'/e) 
performed using the 2001 data sample, recorded 
in somewhat different experimental conditions by 
the NA48 experiment. 

After the 1999 data-taking period, the drift 
chambers of the experiment were damaged by the 
implosion of the beam tube. The data taking 
in 2001 took place with rebuilt drift chambers. 
Thanks to the possibility of a better SPS duty cy- 
cle, the data could be taken at a 30% lower beam 
intensity, allowing the insensitivity of the result 
to intensity-related effects to be checked, and the 
statistics for the final e'/e measurement by NA48 
to be completed. The statistics accumulated dur- 
ing the 2001 data-taking period is roughly half of 
the total statistics accumulated in the 1998 and 
99 periods. The details of the analysis of this data 
set can be found in || . 

1.2. NA48 Method 

The measurement of R proceeds by counting 
the number of events in each of the four de- 
cay modes. The experiment is designed to ex- 
ploit cancellations of systematic effects contribut- 
ing symmetrically to different components of the 
R. Data are collected simultaneously in the four 
decay modes, cancelling the absolute fluxes and 
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minimising the sensitivity of the measurement to 
accidental activity and variations of the detection 
efficiency The Kl and K$ decays are provided 
by two nearly collincar beams. To minimise the 
difference in acceptance due to the large differ- 
ence in average decay lengths, only Kl decays 
occurring in the region also populated by Ks are 
used. The Kl events are furthermore weighted 
flS <X function of the proper lifetime in order to 
equalise the decay vertex distribution of Kl and 
Ks events. With this procedure, the acceptance 
correction cancels at first order in R, thus min- 
imising the systematic uncertainty from Monte 
Carlo modelling of the experiment. 

2. BEAMS AND DETECTOR 

2.1. Beams 

The neutral beams are derived from 400 GeV/c 
protons extracted from the CERN SPS. Because 
of the different mean decay lengths of Kl and 
Ks, two different production targets are used, lo- 
cated 126 m and 6 m upstream of the beginning 
of the decay region. For each SPS pulse (5.2 s 
every 16.8 s), «2.4xl0 12 protons hit the Kl pro- 
duction target. Three stages of collimation are 
used to define the Kl beam, at a production 
angle of 2.4 mrad. Part of the non-interacting 
protons impinge on a bent silicon mono-crystal. 
A small fraction undergoes channelling and pro- 
duces a proton beam of «5x 10 protons per pulse 
transported to the Ks production target, where 
a second neutral beam is derived, at a production 
angle of 4.2 mrad. The Ks beam enters the fidu- 
cial decay region 6.8 cm above the Kl beam. The 
beams converge with an angle of 0.6 mrad and 
the axes of the two beams cross at the position of 
the electromagnetic calorimeter. To distinguish 
Kl and Ks decays, the protons directed to the 
Ks target are detected by an array of scintilla- 
tion counters which comprise the tagging detec- 
tor. The time of each proton is recorded and is 
compared with the time of the decay as measured 
in the main detector. The presence (absence) of a 
proton in coincidence with the event defines the 
event as Ks (Kl)- The Ks beam traverses an 
anti-counter (AKS), formed by a set of scintil- 
lation counters following a 3 mm thick iridium 



crystal. This detector provides an accurate defi- 
nition of the beginning of the decay region where 
it is located by vetoing Ks decays occurring up- 
stream. 

2.2. Main detector 

Charged pion decays are measured by a mag- 
netic spectrometer comprised of four drift cham- 
bers and a dipole magnet giving a momentum 
kick of 265 MeV/c. The momentum resolution is 
a p /p = 0.48% 8 0.009% x p (p in GeV/c). Two 
plastic scintillator hodoscope planes are located 
after the spectrometer. They are used to deter- 
mine the event time of the charged events for the 
tagging procedure. 

A quasi-homogeneous liquid krypton electro- 
magnetic calorimeter with a projective tower 
readout is used to measure the photons from 7r°7r° 
events. The 13212 readout cells have each a cross- 
section of w 2x2 cm 2 . The energy resolution is 
a(E)/E = 0.032/ \J~E~ 0.09/.E 0.0042, where 
E is in GeV. The spatial resolution is better than 
1 mm above 25 GeV. This detector is also used to 
measure the time of the photons for the tagging 
procedure. 

Located after the electromagnetic calorimeter 
are a iron-scintillator hadron calorimeter, fol- 
lowed by a muon counter consisting of three 
planes of scintillators sandwiched between 80 cm 
thick iron walls. 

2.3. Trigger 

The electromagnetic calorimeter is used to trig- 
ger on 7r°7r° events. At the trigger level, the 
calorimeter data are reduced to x and y projec- 
tions which are used to reconstruct the total en- 
ergy as well as to estimate the decay position of 
the event. The efficiency of the 7r°7r° trigger is 
(99.901 ± 0.015)% and is K S -K L symmetric. A 
two level trigger is used for n + Tr~ decays. At 
first level, the hodoscope is placed in coincidence 
with a total energy condition defined using both 
calorimeters. The second level trigger uses infor- 
mation from the drift chambers to perform a fast 
event reconstruction. The efficiency of the 7r + 7r~ 
trigger is (98.697 ± 0.017)%. This is higher by 
0.9% than for the 1998-99 data thanks mostly to 
the lower beam intensity. The efficiency is mea- 
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sured separately for Ks and Kl and the double 
ratio is corrected for the small difference. 

3. EVENT RECONSTRUCTION AND 
SELECTION 

K — > 7r°7r° decays are selected using only data 
from the LKr calorimeter. The details of the re- 
construction and of the cuts applied can be found 
in H and ||. From the measured energies and 
impact point positions on the calorimeter of the 
four showers, the decay vertex position along the 
beam axis is computed assuming that their in- 
variant mass is the kaon mass. The invariant 
mass of the two photons pairs are then computed 
(the resolution is better than 1 MeV/c 2 ) and com- 
pared to the nominal 7r° mass, constructing a x 2 
variable. To reject the residual background from 
Kl — > 37T° events, a cut on x 2 is applied. This 
residual background fraction is (5.6 ± 2.0) x 10 -4 
in the Kl sample, while the Ks sample is back- 
ground free. The background from Ks decays 
produced by scattering of beam particles in the 
collimators in the Kl beam is (8.8 ± 2.0) x 10~ 4 . 

K — ► 7r + 7r~ decays are reconstructed from 
tracks in the spectrometer. In the 7r + 7r~ mode, 
both the longitudinal and the transverse decay 
vertex position can be reconstructed, allowing a 
clean identification of Kg and K L decays. A cut 
is applied on the ratio of the two track momenta 
to remove asymmetric decays in which one of the 
tracks could be close to the beam tube where the 
Monte Carlo modelling is more critical. To reject 
background from semileptonic Kl decays, events 
with tracks consistent with being either an elec- 
tron (from the E/p ratio of the energy deposited 
in the LKr calorimeter over the track momenta) 
or a muon (using hits in the muon counters) are 
rejected. Kinematical cuts are also applied on the 
7r + 7r~ invariant mass (the resolution is typically 
2.5 MeV/c 2 ) and on the reconstructed transverse 
momentum. The residual background in the Kl 
sample is (14.2 ±3.0) x 10 -4 , while it is negligible 
in the Ks sample. 

The fiducial ranges in kaon energy Ek and in 
proper time r used to count events are chosen to 
be 70 < E K < 170 GeV and < r < 3.5 t s , 
where t = is defined at the position of the 



AKS counter and ts is the Ks mean lifetime. 
For Kl events, the decay time cut is applied on 
reconstructed r, while for Ks events the cut at 
r = is applied using the AKS to veto decays 
occurring upstream. The determinations of the 
kaon energy, the decay vertex and the proper time 
in the 7r°7r° mode rely on measurements of the 
photon energies and positions with the calorime- 
ter. The uniformity of the calorimeter response 
is optimised using K e ^ decays and checked using 
7T° and r\ decays produced during special ("??") 
runs in which a 7r~ beam strikes two thin tar- 
gets located near the beginning and the end of 
the fiducial decay region. The absolute energy 
scale is adjusted, with a 0.03% accuracy, using 
Ks — > n°TT° decays, such that the reconstructed 
AKS position matches the true value. It is also 
checked using data from the r\ runs (the r\ mass 
is taken from ||). Taking all uncertainties into 
account (including non-linearities in the energy 
response), the total systematic uncertainty on R 
from the reconstruction of tt tt events is found 
to be ±5.3 x 10~ 4 ||. The uncertainty on R from 
the reconstruction of ir + ir~ events is ±2.8 x 10~ 4 . 

A decay is labelled Ks is a coincidence is found 
within a ± 2 ns interval between its event time 
and a proton time measured with the tagger. Fig- 
ure [l] shows the time distribution for Ks and Kl 
decays to ir + ir~ which have been identified by 
their vertex positions. From this, the misidenti- 
fication probabilities in the 7r + 7r~ mode can be 
deduced: a SL = (1-12 ± 0.03) x 10" 4 for the K s 
tagging inefficiency and ocls — (8.115 ± 0.010)% 
for the Kl mistagging as Ks due to an acciden- 
tal coincidence between the event and a proton. 
As the same tagging procedure is used for 7r°7r° 
and 7r + 7r~ events, R is only sensitive to differ- 
ences in misidentification probabilities between 
the 7r°7r° and ir + n~ modes. Several methods 
have been developed to derive these differences di- 
rectly from the data Q . The results are AasL = 
(0±0.5) x 10~ 4 and Aa LS = (3.4±1.4) x 10~ 4 cor- 
responding respectively to an uncertainty on R of 
±3.0 x 10~ 4 and a correction of (6.9 ±2.8) x 10~ 4 . 
The origin of Aa^s is related to a higher loss of 
7r + 7r~ events from accidental activity in the beam 
and can be predicted using a technique of overlay- 
ing "random" events (taken proportionally to the 
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beam intensity) with tttt decays to estimate these 
losses. The agreement between the observed and 
the predicted Acvls values is illustrated in Fig- 
ure n 
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Figure 2. Measured compared to predicted values 
of Auls as a function of the time during the spill. 
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Figure 1. Time coincidence for Ks and Kl 
7r + 7r _ decays, identified by their reconstructed 
vertex position. 



4. ACCIDENTAL EFFECTS 

The 2001 data were taken at lower beam in- 
tensity to reduce the uncertainties related to ac- 
cidental effects. The overlap of extra particles 
related to kaon decays in the high intensity Kl 
beam with a good event may result in the loss of 
the event. The effect on R is minimised by the 
simultaneous data collection in the four channels. 
The possible residual effect on R can be separated 
into two components: 

1) A difference between the beam intensities 
seen by K$ and K L events: AR = AX x AI/I, 
where A A is the difference between the event 
losses in the ir + ir~ and 7r°7r modes and AI/I the 
difference between the K l beam intensity seen by 
Kl and K$ events. AA is minimised by apply- 
ing to all events the recorded dead time condi- 
tions. The largest one is an "overflow" condition 



in the drift chamber readout, which rejects 11% 
of the 7T7T events (this is significantly smaller than 
the 20% loss in the 1998-99 sample thanks to the 
lower beam instantaneous intensity and a lower 
noise in the drift chambers). A A is estimated 
mainly from the overlay technique and is found to 
be (1.0±0.5)%. AI/I is measured from the rate of 
out-of-time LKr clusters and tracks in good Ks 
and Kl decays. This is illustrated in Figure ||. 
Beam monitors, which have been improved for the 
2001 data-taking period, are also used as cross- 
check. The result is AI/I = (0±1)%. The uncer- 
tainty on R from this effect is thus ±1.1 x 10~ 4 . 
This is significantly better than for the 1998-99 
sample, thanks to the lower intensity and to bet- 
ter beam intensity monitors. 

2) A difference in illumination between Ks and 
Kl decays coupled to a variation of the event loss 
with the impact points of the K° decay products. 
This effect is also computed from the overlay sam- 
ples. No effect is found and the uncertainty on R 
is ±3.0 x 10~ 4 , from the statistics of the overlay 
samples. 
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Figure 3. Probability of accidental activity in the 
LKr in the « 150 ns readout window as a function 
of the time during the spill, for K$ and Kl decays 



5. RESULT AND CONCLUSIONS 

Table [l] shows the summary of the corrections 
and systematic uncertainties on R for 2001 data. 
The residual acceptance correction is related to 
the small angle between the Kl and K$ beams 
and is computed using a large statistics Monte 
Carlo sample. As it is mostly given by the beam 
geometry, it does not rely on a detailed simula- 
tion of the detector. As expected from the design 
of the experiment, all corrections are small. Some 
systematic uncertainties are directly given by the 
statistics of the control samples used to study 
them. From the 2001 sample (comprising 1.55 x 
10 6 K L -> 7r°7r° events, 2.16 x 10 6 K s -► 7r°7r°, 
7.14 x 10 6 K L -> ir+n- and 9.61 x 10 6 K s -> 
TT+TT"), the result R = 0.99181±0.00147±0.00110 
is obtained (where the first error is statistical and 
the second systematic) . The corresponding value 
ofRe(e'/e) is Re(e'/e) = (13. 7± 2.5 ± 1.8) x KT 4 
The agreement between this new and the ear- 
lier results is particularly significant since they 
were obtained from data taken at different av- 
erage beam intensities. Taking into account the 
correlated systematic uncertainty of ±1.4 x 10~ 4 , 



the final combined result from the NA48 experi- 
ment is Re(e'/e) = (14.7 ±2.2) x 10~ 4 . 



Table 1 

Corrections and systematic uncertainties on 



R 



in 10 



=3" 



7T + 7r trigger 


+5.2 


±3.6 


(stat) 


AKS inefficiency 


+1.2 


±0.3 




. of ^°7r° 
Reconstruction r , 

OI 7T 7T 




± 5.3 
± 2.8 




„ , , to 7r°7r° 
Background , , 

tO 7T + 7T 


-5.6 
+14.2 


± 2.0 
± 3.0 




Beam scattering 


-8.8 


±2.0 




Accidental tagging 


+6.9 


±2.8 


(stat) 


Tagging inefficiency 




±3.0 




. statistical 
Acceptance 

systematic 


+21.9 


± 3.5 
± 4.0 


(stat) 


Accidental effect m tensity 

illumination 




±1.1 
±3.0 


(stat) 


Ks in time activity 




±1.0 




Total 


+35.0 


±11.0 
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